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ABSTRACT. 

Held has proposed a coordinate- and gauge-free integration procedure within the GHP 
formalism built around four functionally independent zero-weighted scalars constructed 
from the spin coefficients and the Riemann tensor components. Unfortunately, a spacetime 
with Killing vectors (and hence cyclic coordinates in the metric, and in all quantities 
constructed from the metric) will be unable to supply the full quota of four scalars of this 
type. However, for such a spacetime additional scalars are supplied by the components of 
the Killing vectors; by using these alongside the spin coefficients and the Riemann tensor 
components we have the possibility of constructing the full quota of four functionally 
independent zero-weighted scalars, and of exploiting Held's procedure. 
As an illustration we investigate the vacuum Type N spaces admitting a Killing vector 
and a homothetic Killing vector. We show how the properties of separability, redundancy, 
decoupling and reduction of order can be understood and interpreted in a very general 
manner in our approach, and the advantages of being able to postpone the explicit coor- 
dinate choice to the very last step, when it can be used to simplify the residual ordinary 
differential equations. In a direct manner, we reduce the problem to a pair of ordinary 
differential operator 'master equations', making use of a new zero- weighted GHP operator. 
By first rewriting the master equations as a closed set of complex first order equations, 
we reduce the problem to one real third order operator differential equation for a complex 
function of a real variable — but with still the freedom to choose explicitly our fourth 
coordinate. It is then easy to see there are a whole class of coordinate choices where the 
problem reduces essentially to one real third order differential equation for a real function 
of a real variable. An alternative, more algorithmic approach, using a closed chain of real 
first order equations for real functions, reduces the problem to the same order, but in a 
more natural and much more concise form. It is also outlined how the various other third 
order differential equations, which have been derived previously in work on this problem, 
can be deduced from our master equations. 
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1. Introduction. 

In a previous paper [1] we developed and illustrated a coordinate-dependent integration 
procedure in the GHP formalism [2], but, in spite of its advantages, especially in effi- 
ciency, over the analogous NU [3,4] approach in the NP formalism [5], we feel that it does 
not exploit the GHP formalism to its full potential. We also noted there an alternative 
coordinate-independent, operator approach advocated and developed by Held [6-10], and 
modified and illustrated recently in [11]. In this present paper we will develop this ap- 
proach further — specifically for spacetimes admitting at least one Killing vector — and 
illustrate it by an application to the same problem as considered in the previous paper, [1] 

- the G2 class of the twisting type N vacuum spaces (NT spaces). We emphasise however, 
that both papers are self-contained. 

In the remainder of this section we will outline this GHP operator-integration procedure 

- first of all in general, and then specialised to the particular problem to be considered. 

1.1. Outline of the integration procedure: the general case. 

As discussed in detail in the introduction in [1], the GHP formalism consists of a com- 
plete system of three sets of equations: the GHP commutator equations, the GHP Ricci 
equations, and the GHP Bianchi equations. A procedure for integrating this system was 
summarised in five steps in [11]: 

The first two steps are a coordinate- free integration procedure for the operator p, exactly 
as proposed by Held, [6] ; this p-integration process is a generalisation of the r-integration 
process in the NP formalism. The third step involves the application of the commuta- 
tor equations to three complex quantities — two zero-weighted complex quantities which 
supply four functionally-independent zero-weighted real quantities and one weighted (by 
which we shall mean neither weight being zero, i.e. s 7^ 7^ t or equivalently p 7^ ±g, in the 
usual notations, [2]) complex quantity — so that the commutator equations are replaced 
by an alternative set of equations; we shall refer to these four real functionally-independent 
zero- weighted quantities as 'coordinate candidates' because, at the last two steps, they will 
usually be the obvious (but not necessarily always the most suitable) choice for the four 
coordinates, whereas the one weighted complex quantity will usually be transformed to 
unity by choice of gauge. We can refer to this new set of equations as the 'GHP metric 
candidate equations', but emphasise that although this new set of equations is constructed 
in a manner analogous to the NP metric equations in the NP-NU procedure, and to the 
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GHP metric equations in the earlier paper, [1] there is a fundamental difference; these 
coordinate candidates are chosen from within the system of equations, motivated only by 
considerations of their mathematical structure, whereas in [1] choices were imposed from 
the outside, motivated largely by physical and external geometric considerations. Another 
important difference is that many equations in this new set turn out to be identities, mod- 
ulo the equations in the other sets; in fact we deliberately choose our coordinate candidates 
with the intention of achieving such structural simplifications. 

After these three steps we will have reduced the complete system to a much smaller suf- 
ficient subsystem of differential (operator) equations; this subsystem will consist of (the 
equivalent of) six real tables for the action of the four GHP operators on each of the four 
(real) coordinate candidates and on the one complex weighted quantity, together with any 
residual differential equations from the original system of equations. 

The last two steps involve the introduction of an explicit coordinate system and the adop- 
tion of a specific tetrad gauge; these choices will be made in such a manner that the reduced 
sufficient subsystem of differential equations become as simple and manageable as possible. 
If we choose our four coordinate candidates to be our coordinates, then the four real tables 
for these quantities simply become the definitions of the four differential operators (equiv- 
alently, the tetrad vectors) in this choice of coordinate system; if we choose to use our 
gauge freedom to reduce the weighted complex quantity to unity, then the complex table 
for this quantity simply becomes the definitions of the badly behaved NP spin coefficients. 
Almost inevitably there will be residual differential equations, which we would hope would 
be in a reasonably simple form in the chosen coordinate system; However, this often may 
not be the case, and we may then wish to modify our final coordinate choice to make the 
differential equations more manageable. 

The choice of the four coordinate candidates in step 3 is a crucial, and difficult step. Held 
[6] originally envisaged an 'optimal situation' where (i) six real quantities — more precisely, 
four functionally independent zero- weighted quantities and one complex weighted quantity, 
[11] - are suggested explicitly by the spin coefficients and Riemann tensor components, 
and (ii) the integration procedure, after the first three steps, yields only a complete invo- 
lutive set of tables of the GHP operators on these six quantities. (In such a situation the 
problem would be essentially solved; since the differential operators — equivalently the 
tetrad components — could be written down directly from the tables with the four zero- 
weighted quantities as coordinates, and the complex weighted quantity gauged to unity.) 



4 



Unfortunately, so far, in any practical application of this method, it has not been possible 
to obtain, directly from the spin coefficients and Riemann tetrad components, the required 
four functionally independent zero- weighted quantities. In practice, less than the full quota 
of four zero-weighted quantities is usually supplied by the spin coefficients and Riemann 
tensor; at this stage it has been customary [12] to make a translation back to the NP 
formalism or to introduce the 'missing coordinates' from outside the formalism [6,7,10,13]. 
So, in fact, usually a compromise has had to be made at this stage, and the ideals of the 
coordinate-independent integration procedure in the GHP formalism as proposed by Held 
have had to be modified. However, we shall show in this paper that it is possible, even 
when the spin coefficients and Riemann tensor components fail to provide the full quota 
of zero-weighted quantities, to continue working in the GHP formalism along the lines 
proposed by Held. 

We now know why, in practice, it has not been possible to obtain the full quota of four 
functionally independent zero-weighted quantities from the spin coefficients and Riemann 
tensor components: the investigations in GHP formalism which have been carried out were 
in specialised classes of spaces, which almost inevitably means that there exists at least 
one Killing vector. This in turn means the existence of at least one cyclic coordinate in 
the metric, and hence our inability to obtain four functionally independent quantities from 
quantities constructed from this metric (e.g. from spin coefficients and Riemann tensor 
components). However, recently, without going outside the GHP formalism, a means of 
overcoming this problem has been proposed and illustrated with an example in [11]. The 
space under consideration (implicitly) contained two Killing vectors, and it was found 
that the spin coefficients and Riemann tensor components only supplied two functionally 
independent scalars; the additional two coordinate candidates were obtained by taking 
zero-weighted 'potentials' for some of the spin coefficients and Riemann tensor components. 
The procedure for picking out such a potential relies on intuition rather than any standard 
procedure; we would prefer a more algorithmic method where various possible choices of 
coordinate candidates can be generated and tested directly, and the most suitable chosen. 

In this paper we present such an alternative method for choosing coordinate candidates. 
Fortunately, in spaces containing Killing vectors, there is another source from which we may 
find additional coordinate candidates; if we include explicitly the Killing vector equations 
alongside the other field equations, then — since the Killing vector components are not 
constructed from the metric — we can look to the Killing vector components for additional 
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scalars, functionally independent of the metric. So we have found another possible way of 
satisfying the first condition of field's optimal situation — although to obtain the required 
six quantities we have to go beyond the spin coefficients and Riemann tensor components, 
which was field's original proposal. This is precisely the situation which occurs in the 
G2 class of the NT spaces, and makes it an ideal application to illustrate this method. 
(We should emphasise that we are not claiming, in this paper, that the Killing vector 
components will always supply all the missing coordinate candidates; the relationship of 
zero-weighted GHP quantities to Killing and homothetic Killing vectors needs careful and 
detailed treatment, which will be given elsewhere, in this paper, we are simply pointing out 
the possibility of using these components, and giving an application where it is possible.) 

ft is emphasised that once the tables are found for the derivatives of our four coordinate 
candidates and for the one complex weighted quantity, then — alongside the Ricci, Bianchi 
and Killing equations — we have all the information from all the equations. We would hope 
that this complete system would reduce to satisfy the second condition of Held's 'optimal 
situation' i.e. that only a complete involutive set of six tables remain. However, this is 
unlikely in practice; we can expect that in the tables some additional unknown functions 
will occur, which have to satisfy some residual differential equations from the Ricci, Bianchi 
and Killing sets, fn such a situation the choice of the coordinate candidates as coordinates 
is not necessarily the best choice; an alternative choice may make the residual differential 
equations more manageable. At this stage we could of course write down the differential 
equations with the coordinate candidates as coordinates, and attempt to obtain further 
simplifications by explicit coordinate transformations, in the usual manner. However, we 
believe that there are still advantages to be gained by keeping within the GHP formalism, 
obtaining further simplification within that formalism, and finally allowing the structure 
of the suitably simplified GHP operator equations to suggest the optimum coordinate 
choices. Again, the application in this paper illustrates precisely how this procedure can 
be implemented. 

1.2. Outline of the integration procedure: the G2 case in the NT problem. 

For the particular problem being considered in this paper the complete GHP system con- 
sisting of the three sets of equations — Bianchi, Ricci and commutator equations — has 
to be supplemented by three other sets of GHP equations: the Killing vector equations, 
the homothetic Killing vector equations, and the non-Abelian G2 condition. The five steps 
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in the integration procedure will then be carried out for the whole system of six sets of 
equations. Before beginning, we will of course ensure that all the equations are specialised 
to vacuum Petrov Type N spaces. 

Held, in [8,9], has investigated some algebraically special vacuum spaces admitting Killing 
vectors, and in Section 2 and in the Appendix we clarify some differences in approach and 
notation between his work and our previous paper, [1]; in Section 3 we make use of the 
first part of Held's work in [8] - the p— integration — and specialise the relevant calcu- 
lations to Petrov type N spaces, as well as extending to the G 2 case. This p— integration 
(corresponding to the first two steps of the integration procedure) is summarised at the 
beginning of Section 3, and we then also carry out the third step — applying the com- 
mutators to six appropriate quantities. The choice of these quantities is dictated by our 
wish to keep the calculations as simple and manageable, and to end up with as concise a 
system, as possible; however, we also try to keep in touch with the earlier method in [1], 
and to draw comparisons. Since the spaces under consideration have one Killing vector we 
will not be able to obtain all our coordinate candidates directly from the spin coefficients 
and Riemann tensor components; in fact none of our candidates is chosen in this manner: 
three candidates are constructed from the tetrad components of the Killing vector and 
homothetic Killing vector combined with the spin coefficients, while the fourth candidate 
is chosen as a potential for a combination of some of these quantities. Since we choose 
our six quantities as far as possible to suit the structure of the system of equations, and 
to give as simple a presentation as possible, when the tables are constructed for these six 
quantities we find considerable redundancy and simplifications once these equations are 
put alongside the Ricci, Bianchi and Killing equations. Unfortunately, we do not achieve 
Held's optimal situation, since in addition to the six quantities, the six tables involve other 
functions, which are themselves subject to two residual Ricci equations. 

The fourth and fifth steps of the integration procedure involve the choice of coordinates and 
gauge. In the last part of Section 3 we make the most straightforward choice — choosing 
the coordinate candidates as coordinates. We find that these choices do cause considerable 
structural simplification; in particular they enable the residual partial differential equations 
to separate into a pair of coupled ordinary differential equations for one unknown complex 
function of one real coordinate. However, these residual ordinary differential equations do 
not look very manageable, and, in particular, decoupling seems a problem. 

In our work in Section 3 our primary concern is to ensure that we consider a sufficient 
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system of equations; up to this stage we are not so concerned with the most efficient 
presentation of the residual partial differentiatial equations. However, we do keep in mind 
the need for subsequent simplification, and of course we are looking for the simplest options; 
and in particular we are hoping for candidates which permit separation. So although, in 
this case, we will ultimately choose coordinates different from the coordinate candidates, 
yet we realise that they will have to be closely related in order to share the separation 
property; specifically we will retain the first three coordinate candidates as coordinates, but 
we will allow the structure of the residual pair of ordinary operator differential equations 
to suggest our fourth coordinate. 

The great advantage of our GHP approach is that we can get an overall picture of the 
structures of the equations, and this leads to insights not just into separability but also 
into redundancy, decoupling and order reduction in a manner which does not first require a 
precise coordinate choice; the final coordinate choice can then be made in an informed man- 
ner which best exploits the separability and decoupling properties, and seeks for maximum 
order reduction and/or simplification in the final differential equation. These advantages 
are demonstrated and illustrated in the remaining sections. 

In Section 4 we go back to the residual pair of Ricci equations as given in the GHP for- 
malism, and by modifying the GHP weighted operator f> to a new zero- weighted operator 
V we obtain a simpler version of these operator equations, where it is explicit that all the 
terms depend on only one real variable. We consider this version of these equations as our 
'master equations'. From these equations we are able to experiment with various different 
approaches to decoupling and reduction of order; in particular, we find one particularly 
structurally simple presentation of the master equations, as a set of first order ordinary 
differential operator equations in complex functions. 

In Section 5, we show that by choosing our fourth coordinate from a class which exploits 
the structure of this simple version of the master equations, the decoupling problem is less 
complicated than in the previous coordinate choice in Section 3; by a little manipulation the 
problem reduces to one real third order differential operator equation for a complex function 
of a real variable. However, we still have the freedom to make our explicit coordinate choice, 
and by a suitable choice the problem easily reduces to one real third order differential 
equation for a real function of a real variable. We note that there is a class of such 
coordinates which reduce the problem to this order. 

In Section 6 we present an alternative and more algorithmic approach to obtaining the 
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residual equation. By a systematic approach we construct, from our set of complex first 
order equations, a set of real first order equations in the form of a closed chain. These equa- 
tions suggest both the cordinate and dependent variable in a natural way; the decoupling 
problem then disappears, and we obtain a comparatively concise form for the residual real 
third order differential equation. 

The early work on the G 2 case of the NT problem by Mcintosh [14] reduced the problem 
to a complicated third order complex differential equation for a complex function of a real 
variable, which he rearranged to a sixth order complex differential equation for the complex 
function alone — without its complex conjugate occuring explicitly. In recent years the 
problem has been reduced to essentially a single real ordinary differential equation of a real 
function of a real variable; the order of the final differential equation has gradually been 
lowered (at the expense of increasing non-linearity) and most recently there have been 
obtained a few (very complicated non-linear) third order real differential equations for a 
real unknown function, each highly dependent on its background formalism and particular 
coordinate choice [15-18]. The only known solution, the Hauser solution [19,20,21], is found 
to be a singular case. Each of these analyses began with a preferred coordinate system, and 
a lot of involved manipulation and complicated coordinate changes were needed in order 
to arrive at the final third order equation; it is not clear from each individual case whether 
alternative (and perhaps simpler) third (or lower) order equations could be obtained, nor 
is there any obvious links between the very different equations so far obtained. On the 
other hand, our method supplies us with a whole class of third order equations, and with 
the possibility of obtaining even more in a manner where we have some understanding of 
the structures involved, as well as some freedom to simplify these structures. 
We highlight the insights we have obtained into the GHP integration procedure and into 
this particular application in the concluding section; we also outline how the various dif- 
ferential equations obtained previously by other approaches to this application can be 
deduced from our master equations. 

2. Formalism and Notation. 

As noted above, Held [8,9] has used the GHP formalism to investigate algebraically special 
vacuum spacetimes admitting Killing vectors, and we shall make direct use of some of his 
results in this work. Some points need to be clarified at the outset: 
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Firstly, Held, [6,8,9] slightly modifies three of the usual GHP operators p', <9, & to f> , d, d ; 
this makes no essential difference to the overall structures of the GHP formalism, but sim- 
plifies calculations by taking advantage of some properties of algebraically special vacuum 
metrics. In particular, it is easy to carry out a 'coordinate- free p-integration'. 

Secondly, as noted in Section 1, the set of Killing equations with which Held [8,9] works 
differ a little — being somewhat simpler — from the conventional set, and in particular 
from the set used in [12], and in the previous paper. We show in Appendix I the relationship 
between these two sets, and confirm that the two apparently different sets of equations in 
the two papers are equivalent for Killing vectors. We also generalise the set in [8,9] for the 
presence of homothetic Killing vectors, and show its equivalence to the analogous set in 
[12]. 

Thirdly, there is a question of notations, and this needs to be set out in detail to avoid 
confusion: 

(i) In the previous paper, [1] we used the symbols 6, c, a for tetrad components of Killing 
vectors, and E° for the twist. Held uses the symbols £ 0) £i) £2, [9] or 6> , #i, 6> 2 , [8] for Killing 
vector components, and 0°(= — 2i£°) for the twist. When we make use of Held's results 
we shall translate them into the former symbols for consistency with the previous paper. 

(ii) Quantities which are annihilated by the operator p are labelled in Held's work with 
e.g. p?] = 0; an equivalent notation has already been adopted in the previous paper. 

(iii) Standard GHP usage, also followed by Held, uses the prime notation for half of the spin 
coefficients, and half of the GHP operators; in this paper we shall also use that notation. 
Although our previous paper uses the more familiar NP versions of the spin coefficients, 
since only k'(= —v) and p'(= —/j,) occur explicitly in this paper there will be no difficulty 
in comparison. 

(iv) The letter P has a special usage in NP notation [3,4], and has been used extensively 
in the previous paper, [1]; the related P° plays an important role in Held's work [8,9,10], 
and we retain it in this paper. If a comparison is being made we must take into account 
a factor of 2, (P° ~ 2P) as can be seen in the defining equations in each paper. We also 
point out that P defined in [1] is real. 
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with 



3. An operator-integration approach with coordinate candidates as coordi- 
nates. 

3.1. Steps 1,2: The p— integration. 

We specialise the p— integration results in [8] to the Petrov type N case with the substitu- 
tions 

^2=0 = ^3 (3.1) 

giving: 

The residual Bianchi and Ricci equations. 

p' = pp'° (3.2a) 
k' = k'° (3.26) 
* 4 = pn (3.3) 

= _dV° (3.4a) 

& p'° = 2zEV° (3.46) 

cV° = £ V ° (3.4c) 

dd'z° = 2E°p° (3 Ad) 

PP = P 2 i>'p = P 2 P'° 

dp = dp = -2ip 2 dZ° (3.5) 

The commutator equations. 

[p, ~d] = (3.6a) 

[p,d'] = (3.66) 

[£,£']= (3.6c) 

[i>\d] = -— $> + qR'° (3.6a 1 ) 
P 
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and 



&',&'] = -y$>+ P K'° (3.6e) 

[d, ~d'] = {P-- - 2zE°f' + pp'° - qp'° (3.6/) 

P P 

(In addition to the Bianchi and Ricci equations given above there are the two equations 

2z£ E° = + p'° 

= o 

which follow from [8]. However, both are identically satisfied by virtue of the other Bianchi 
and Ricci equations, given above, and the commmutators. We show this as follows: apply 
commutator (3.6f) to E° and use (3.4d) and its conjugate to get the first; apply commutator 
(3.6f) to k'°, followed by (3.6e) applied to p'° and then use appropriate equations from 
(3.4) together with the first equation to get the second.) 

The equations for one Killing vector. 
When the Killing vector £i is given by 

= ai l» + bm* - Cl m^ - c x m» (3.7) 

the Killing equations can be p-integrated to obtain, [8] 

a x = -lic^Y? - p' ^ + zE°£ 'b° + + ht>'b° 

2 p p 

h = bl ( 3 - 8 ) 

a = c\lp 

with 

£'&? = -i + (3.9a) 

~db\ = -2iE°c? (3.96) 

f 'c\ = = del (3-9c) 

»ipV° ~ <$p'° - c\dp'° + 2p'°^b\ = (3.9d) 
b°(p'° ~ P°) + 2ic?aE° + 2ic p 1 d'x° + zE°(3c? + d'cl) = (3.9e) 
where the subscript i denotes quantities associated with the Killing vector £i. 
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The last equation was not explicitly displayed in [8] but is given in [9] where it is pointed 
out that it follows from the condition that a\ is real. 

However, when we examine this last equation we note that it is simply the commutator 
(3.6f) applied to b°. In addition, when we apply the commutator (3.6d) to b° and differen- 
tiate the resulting equation by d , after using (3.6f) a number of times together with some 
equations from (3.4), we obtain (3.9d). Hence the two equations (3.9d,e) are identities 
modulo the residual equations and the commutators, and so can be omitted. 

The equations for one homothetic Killing vector. 
When the homothetic Killing vector £2 is given by 

f£ = a 2 l tI + b 2 n tI - c 2 m^ - c 2 m^ (3.10) 

we use the set of equations (AI.2) in the first Appendix, and in the same manner as in [8] 
we can easily integrate to obtain 

a 2 = -2ic° 2 ~d'z° - p'°b° + iE°i>'b° 2 + + -)(i>'b° 2 - </>) 

2 P P 

b 2 = b° 2 (3-H) 

C2 = C° 2 /p 

with 

^b 2 [ = -\{dc° 2 + ~d'c° 2 - ( j ) ) (3.12a) 

~db° 2 = -2iY,°c° 2 (3.126) 

f>' c ° = = d'c° 2 (3.12c) 

b°$>'p'° - c° 2 d' P '° - c°dp'° + 2p' ^b° 2 = (3.12d) 

b 2 (P° ~ P°) + 2ic0Z° + 2ic°d'z° + iZ°(dc° 2 + ~d'c° 2 - 0) = (3.12e) 

where the subscript 2 denotes quantities associated with the homothetic Killing vector 
£2- We note that in this case also, the last two equations (3.12d,e) are identities modulo 
the residual equations and the commutators, and so can be omitted. 

The non-Abelian G 2 condition. 
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We now consider the case for a non-Abelian G 2 of homothetic motions; by an appropriate 
choice of basis in the Lie algebra we get 



[£i,£ 2 ] + 2^ = 



(3.13) 



where £i(= £fV M ) and £ 2 (= ^2^m) are the Lie derivative operators associated with the 
Killing vector £1 and the homothetic Killing vector £2 respectively. Equation (3.13) has 
to be satisfied when applied to an arbitrary (scalar) quantity or equivalently to four real 
functionally independent quantities. In the present context that means that it must be 
applied to four real zero-weighted functionally independent quantities, and for such zero- 
weighted scalars the operators can be written as 



So for the case of Petrov type N vacuum spaces admitting one ordinary Killing vector £1, 
and one homothetic Killing vector £ 2 , we would be required to solve the residual Ricci 
and Bianchi equations (3.4,5) together with the residual Killing equations (3.9a,b,c), the 
residual homothetic Killing equations (3.12a,b,c), the commutator equations (3.6) and 
non-Abelian condition (3.13). 

3.2. Step 3: Choosing the coordinate candidates and applying the commuta- 
tors. 

Since the first two steps of the procedure set out in [11] have been carried out, the next 
step is to ensure that the commutator equations are completely satisfied; to ensure this we 
must apply them to four real, zero- weighted, functionally independent quantities and to one 
complex weighted quantity. We have just noted the similar requirements for (3.13) so it will 
be convenient to apply both sets of equations — (3.6) and (3.13) respectively — to the same 
four coordinate candidates. We would prefer to select those four coordinate candidates in 
a manner which would ensure that any additional constraint equations resulting from the 
commutator equations assume as simple a form as possible; but we also wish at this stage 
not to move too far from the coordinate choices of the first method. Fortunately, for the 
first three coordinate candidates, these two wishes coincide. 

For each of the four real and one complex quantities we will obtain a table for the action of 
the operators; we may also obtain some new differential equations when the commutators 



i x {= £! V,) = a lP + 
£ 2 (= ^V l ) = a 2 p + 6°f 



c° 2 ~d - c° 2 ~d' 



(3.14) 
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and G 2 condition are used on each of these quantities, as well as some simplifications. We 
now set out these results systematically for each quantity in turn. 



Choice of the first two real (one complex) coordinate candidates, (. 

Noting the simple equations for the (0,2) quantities c° and c% we begin with the compar- 
atively obvious choice for our first (complex zero- weighted) coordinate candidates 

C = ii^K) (3-15) 

so that 

K = K = ~d'C = (3.16) 

(We choose the factor i in our definition of £ in order to make easy direct comparison with 
[1], and also because this choice simplifies a later part of the calculation. We also omit ° 
on ( , even though it is annihilated by |>. ) 

When the commutators are applied to this zero-weighted quantity the only non-trivial 
results are, 

= i>'p° = dP° = (3.17) 

where P° is the (0,-2) quantity, 

P° = d( (3.18) 
Noting the similarity between these equations and those for c^, we can write 

J>( C °P°) = £ '( C ?P°) = d'(c°P°) = (3.19) 

Since (c°P°) is also a zero-weighted quantity comparison between its equations (3.19) and 
those for ( show that 

c\P° = 2g(Q (3.20) 

where g is an arbitrary function of (. (The factor 2 has been introduced to give correspon- 
dence with [1]) 

Turning next to equation (3.13), and applying (3.14) to U we find 

*iC = -4k = -2.(0 

= -4k = iclCdC = 2*07(0 

and so obtain from (3.13), 

9(0 = i (3-22) 
15 



or 




and 




(3.23) 



where ( satisfies (3.16) and P° satisfies (3.17). 

In summary, we have a table (3.16,18) for the action of the operators on the two real 
coordinate candidates £; and the action of the commutators on these quantities yields the 
additional equations (3.17) — the table giving the action of most of the operators on the 
weighted quantity P°. The result of applying the condition (3.13) has been to obtain 
explicit expressions for the Killing vector components c°, c%- 

(Although this is not immediately apparent from the respective definitions, a little work 
confirms that ( as defined above essentially agrees with £ as defined in [1] . We could have 
chosen to define ( in this paper in a manner more closely analogous to [1] by introducing 
it as a potential for P°, which, in turn, could have been introduced as a potential for some 
of the spin coefficients. However, as we have emphasised before, we wish in this paper to 
illustrate the direct method of choosing coordinate candidates.) 

Choice of the third coordinate candidate, R. 
A rearrangement of (3.5) 





2/\P° 



{ P '° + P'°)/\P 



.0 



iO 



(3.24) 



suggests an obvious choice for our third coordinate candidate R, 



R=-( - + -)/\P 
P P 



(3.25) 




$>'r = (p' + p' )/\p° 




(3.26) 



d'R = -2i 




P' 



o 



-R 



g'po 



2P° 
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From their respective differential equations given in tables (3.16) and (3.26), the three 
coordinate candidates (£ + C)>*(C ~~ C)>-^ are easily seen to be functionally independent 
of each other; and when the commutators are applied to R, the only two commutators 
(3.6d,f) which are not identically satisfied, give respectively, 

*'° = -f'(^) ( 3 - 27 «) 
-/ dP° 

P'° = -d (f^) (3.276) 
Substituting these values into (3.4a) shows that 

~/~/ d'P° 

*4 = <??( — ) (3.28) 
and also that (3.4c) is identically satisfied by virtue of (3.6d). 

So the only remaining Ricci and Bianchi equations to be solved are (3.4b) and (3.4d) which 
now become 

~/~ d'P° ~i d'P° 

v Do ' r v Do ' 

F F 3.29 

~~/ ~/ dp° 

dd e° = -e°<9 - — ) 

v po ' 

The first of these may be rearranged to 



.-./. d'p° 



dd (^— ) = (3.30) 

po 6 l z 

There are also the inequalities 

E V (3.31a) 
~~/ dP° 

9 P (^) * (3.316) 

The remaining constraint equations from the Killing equations and from the homothetic 
Killing equations, after the substitutions (3.23), are 

-/ 8P° d'P° 

r 1 V po 2 po 2 ' 



~db\ = -4E°/P° 
~Bb° = -4iCE°/P° 
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(3.32) 



Turning next to equation (3.13), and applying (3.14) to R we find 

6 (3.33) 

£ 2 i? = (2 + 

and so we easily confirm that (3.13) is satisfied identically when applied to R. 

In summary, we have a table (3.26) for the action of the operators on the coordinate 
candidate R; and the action of the commutators on R yields explicit expressions for some 
of the spin coefficients and ^4; some of the Killing equations are also simplified. The result 
of applying the condition (3.13) to R yields no new information. 

Choice of the fourth coordinate candidate, u. 

We have noted in the introduction to this paper that we can introduce each of the coor- 
dinate candidates either directly as one of the (combination of) existing elements in the 
GHP formalism, or as a 'potential' for existing elements. Such potentials are chosen by a 
careful consideration of the structure of the equations, especially the commutators, pay- 
ing particular attention to the appropriate weights, [11]; specifically the four derivatives 
p?7, p rj, drj, d n of a zero-weighted potential r\ are equated to (combinations of) some of 
the spin coefficients and Riemann tensor components in such a way that the commutators 
acting on r\ are identically satisfied. Since for zero-weighted 77, 

V^ry = n^rj + Z M f> n - rh^drj - m^d'n (3.34) 

clearly fixing its GHP derivatives determines r\ uniquely (up to an addititive constant). 

We can look for hints for possible potentials by rearranging the residual equations, paying 
particular attention to weights; for instance the first and last of (3.32) can be combined to 
give 

~ / - - dP° 4 - 6 

P(^-^C) = (C + C)^2 -~2^ ( 3 - 35 ) 



i.e. 



■P°|P°|(&S-t&fCh _oS|p ,_ (4-0)P o |P°| 



p \ (t + n )~° 



(C + C) > ' ' 2(C + C) 

2(C + C) 1+m a(|P°|(C + C)- 1 - m ) with m=-^ 



(3.36) 



Therefore, the choice 

%>u = (3.37a) 
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i>'u=\P°\(( + 0- 1 - m (3.376) 

~du = P°\P°\(b° 2 - tbIOiC + C)- 2 " m /2 (3.37c) 

where u is a real zero-weighted quantity, guarantees that commutator (3.6a) is satisfied 
identically when applied to u. (In fact it follows easily that (3.6f) is the only commutator 
not satisfied identically when applied to u.) 

We now need to check explicitly that the four coordinate candidates R, £, £, u, are function- 
ally independent. We can write the four derivatives of each of these four real zero- weighted 
quantities in determinant form as 

2/\P°\ (p'° + p'°)/\P°\ (2idZ°/\P°\) -RdP°/2P° (-2id'E° /\P°\) - Rd'p°/2P° 
P° 

P° 

o ip^kc + c)- 1 -^ p^c + cr 1 p^C + C)" 1 



(3.38) 



where we have defined the function (, u) by, 



= (C + Qdu (3 3Q) 

po 



Hence, using (3.39), (3.6d) and (3.37b), we have 

hP° hf 

^2 = 2 j^(C + CT + 2(1 + m)(C + C)" 1 (3-40) 

It is easy to confirm that this determinant is non-zero. (It can be seen that this definition 
for I agrees with the definition in the previous paper, [1].) 

When we apply the remaining commutator (3.6f) to u we obtain 

2|P°|E° = i0u-d'du)(( + () 1+m (3.41) 

which becomes 

2|P°|E° = i(p°M-p°d'e+ |P°| 2 (C + C)"V-£))(C + C) m (3-42) 

We can find 6°, b% explicitly from (3.37c), 

K = 2,(^^ + 0"' (3.43) 
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m 




~p° 


\p°\ 


o'e 
~p° 


\p°\ 



bi = 2 (« + <3(C + <T (3.44, 
and when we substitute these values into the remaining Killing equations (3.32) we get 

-(C + C) m (3.45a) 

-(C + C) m (3-456) 
We can now use (3.39) and (3.37b) to rewrite (3.37c) as 

^ = ^(00- (3-46) 

and a comparison with (3.45) shows that these equations give the very crucial simplification 

that I and u are functionally dependent. 

When (3.45) are substituted into (3.42) we obtain 

2 W\ = 1 ( J^T ~ W\ + (c + c) " 1_m(£ ~ (c + c)2m (3 ' 47) 

The two residual Ricci and Bianchi equations (3.29) now become equations for t by the 
substitutions (3.40,47). We now seem to have too many equations — one real and one 
complex — for the complex unknown t. However, we shall see that this is compensated 
for by the fact that these two equations are not independent; it will be easier to show this 
explicitly when we introduce our coordinates. 

The inequalities (3.31) can also be rearranged with these substitutions. 
Turning next to equation (3.13), and noting that 

£i« = 

(3.48) 

£ 2 u = 

we confirm that (3.13) is satisfied identically when applied to u. 

In summary, we have a table (3.37) for the action of the operators on the coordinate 
candidate u; and the action of the commutators on u gives explicit expressions in terms 
of £ for the twist E° and for the Killing vector components b°, b% : as well as the action 
of the remaining operator d on weighted P°. The remaining two Killing equations give 
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the important result that the complex function £ and the coordinate candidate u are 
functionally dependent. The result of applying the condition (3.13) to u yields no new 
information. 

Choice of one complex weighted quantity, P° . 

When we apply the commutators to P° and all the earlier equations are taken into account, 
all the commutator equations acting on P° are identically satisfied. 

Summary of results in this subsection. 

Having confirmed that our four coordinate candidates are functionally independent, we 
now know that we have extracted all the information available from the commutators and 
the non-Abelian condition. The net result is that we have obtained explicit expressions for 
all the spin coefficients, Weyl tensor component and Killing vectors components, as well as 
six tables of twenty-four (real) equations — (3.16,18) and their complex conjugates, (3.24), 
(3.37a,b,39), (3.17,40) — for the action of the operators on the four coordinate candidates 
.R, u, C, C an d on the one complex weighted quantity, P°. However, this set of tables is not 
self-contained; it contains explicit expressions in £, which is a function of u. The quantity 
I must satisfy the two residual Ricci equations (3.29), subject to the inequalities (3.31), 
when the substitutions (3.40,47) are made. 

The next obvious step is to use the coordinate candidates as coordinates. 

3.3. Steps 4,5: Choice of coordinate candidates R, £, £, u as coordinates and 
choice of gauge. 

We will now choose the coordinate candidates as coordinates R, £, £, it, which means that 
their respective four tables of equations (3.16,18), (3.24) and (3.37) become identities, 
simply defining, in these coordinates, the GHP operators when acting on zero-weighted 
quantities as, 

P = HR) 9r + p(C) d c + P (C) 0f + £(«) d u = 2\P°\~ 1 d R 
p' =(p'° + p'^lPT'dR + \p°\(C + Q- l ~ m d u 

/ r)Y° P°h f \ 

d =(2*^| - ^^j-(c + cr - rp°(i + m)(c + c)- 1 )^ + p°d c + p°£(c + cr'du 
~&' =(- 2 ^ - R ^f^ + o m - rp°(i +m)(c+ cr^dn + p°d- c + p°i(c + o~ i d u 

(3.49) 
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The Killing operator becomes 

£ t = -2id c + 2id^ (3.50) 
while the homothetic Killing operator becomes 

£ 2 = -2(d c - 2C<9 C - - AmRd R (3.51) 

There is of course also the table of operator equations for the weighted P°, but this simply 
yields the badly behaved spin coefficients (if we want them) when we make our choice of 
gauge. 

This coordinate choice exploits the separability property, and we find that the problem has 
been reduced to only two differential equations which can now be written out explicitly in 
these coordinates. Since £ is a function of u only, i.e £ = £(u), we can replace £> £ with 

£l = i£ u = i\po^ +(:) -l-m (352) 

where ' denotes differentiation with respect to u and we have used (3.37). So the only 
remaining equations to be solved are the two residual Ricci and Bianchi equations (3.29) for 
£ which is a complex function of the real coordinate u. When the substitutions (3.42,47) 
together with (3.58) are made into the equations (3.29) and the inequalities (3.31) we 
obtain 

££ 7 + Mil + l(2m + 1)£ - 2££ - 2£ 2 - (Am + 2)1 - 2m[m + 1) = (3.53) 

{£P £ + U{[2m + 1)1+ (2m - S)£ + U£ + ££} + IH + (2m - 3)U 2 
+ (2m + l)££l+ (Am - 2)11 - (4m 2 -Am- 3)££~ - (Am 2 -Am- 2)1} 
- {c.c.} 

= 

££-££ + £-£^0 (3.55) 

£ £ + U£ + (2m + 1)£ ^ (3.56) 

(These correspond to the equations and inequalities (111,112) in the previous paper [1], 
as well as to the equations and inequalities (25,26,28,29) originally given in [17], when the 
latter two sets of equations are specialised to / = constant.) 
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We can substitute (3.53) (and its complex conjugate) into (3.54) to obtain a simpler version, 

{££(! + 2m - 1)£ + ££ 2 I + (2m - l)££ 2 

+ (2m + 1)£££ + (4m 2 - 1)11 - (Am - 2)££ - (6m 2 -2m- 2)£} 

(3.57) 

- {cc} 
= 

We now find explicitly the redundancy relation between this pair of equations to be 

4- [eq. (3.57)1 = (2m - 1 + 1) [eq. (3.53)1 - (2m - 1 + £) [eq.(3.53)l (3.58) 

The gauge choice does not affect the two residual differential equations, but putting P° = 1 
will cause some minor simplification in the form of the differential operators. 
So we have succeeded in reproducing the results in the previous paper [1] by a method 
which was clearly influenced by the choice of the fourth coordinate u in that paper. 

3.4 Summary of this section. 

In Sections 3.1 and 3.2 we have succeeded in reducing the problem to essentially a pair of 
coupled ordinary differential equations for £ by a method which was coordinate invariant. 
It may be argued that we have really introduced coordinates (the coordinate candidates) 
in all but name, so perhaps we should emphasise that our method is coordinate invariant 
in the sense that there is no background coordinate metric imposed, and that we do not 
have to adopt the coordinate candidates as our final coordinates. We emphasise again 
that it is structurally imperative that we apply the commutators to (the equivalent of) 
four functionally independent real scalars to ensure that we have a complete system of 
equations. Therefore, our prime concern was simply to establish a complete system and 
we have not been thinking of the suitability of these quantities in their optional additional 
role as coordinates; although of course we are always concerned with getting the complete 
system of equations in a reasonably concise and manageable form. Hence whether our 
coordinate candidates are the coordinates which give separation, decoupling and reduced 
order is of secondary importance, at this stage; it is of course preferable if it happens, as 
it makes subsequent work shorter. 

Having obtained explicitly a complete system we then decided to tentatively adopt the 
coordinate candidates as our coordinates and to explore their usefulness. We have seen in 
Section 3.3 that this particular coordinate choice does give separation, and this is a bonus 
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for us; of course, we also want the equations to decouple and reduce to as low an order as 
possible. Although we could continue to work in the coordinates chosen in the last section 
we have pointed out before that we have the freedom to choose coordinates other than 
the coordinate candidates; therefore, we will now investigate these possibilities. From now 
on, having learned something of the structure of the residual differential equations still 
to be solved, motivation for preferring a particular choice will be its usefulness not just 
in separation, but also in enabling us to decouple the residual equations and to reduce the 
order of the equations. 

4. The master equations, and alternative coordinates to coordinate candidates 

Although allowing the coordinate candidates to become the coordinates, as in the last 
section, is the most obvious choice, it is not the only one. Of course, in practice, we 
want to choose as coordinates those four quantities on which each of the four derivative 
operators yield simple expressions, since these expressions give us the explicit form of the 
operators and tetrad components in that coordinate system; as well we recognise that 
we will be left at the end of the operator-integration procedure with a residual set of 
differential equations still to be solved, and further progress towards their solution will 
depend on a suitable choice of coordinates. However, the particular choice made in the 
last section enabled us to obtain separation of the coordinates and reduce the problem 
to ordinary differential equations in the fourth coordinate u. Clearly we would wish any 
other coordinate choice to give us the same separation properties; so we will retain the 
first three coordinate candidates as coordinates, but postpone the explicit choice of the 
fourth coordinate. To make a skillful and informed choice of the fourth coordinate we need 
more understanding of the precise structure of the residual differential equations. In this 
section we will first get a better picture of this structure, which will motivate the eventual 
explicit choice for the fourth coordinate; so we will now continue on directly from the end 
of subsection 3.2 — where we had reduced the problem to solving the residual pair of Ricci 
equations (3.29), subject to the inequalities (3.31) — and develop further the equations 
(3.29) in the GHP operator notation. 

4.1. The master equations and redundancy. 

We have noted in Section 3.2 the crucial result that £ and u are functionally dependent. Fur- 
ther, when we compare the differential equations (3.45) for £ with those for u, (3. 37b, 3. 46) 
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we note that 

= V(£)V l u (4.1) 

where 

P(£) = (C + C) 1+m ^ (4.2) 
so £ and T>(£) are functionally dependent on u and hence on each other. The operator V 

(( + C) 1+m ~> 

V = ,po| ( 4 - 3 ) 

is simply a zero-weighted operator formed by scaling the (l,l)-weighted operator k ; this 
new operator also has the property that when it operates on a zero-weighted function of 
u, it yields another zero-weighted function of u. 

We now wish to write out explicitly the two residual Ricci equations (3.29) in this notation, 
and so we first rewrite (3.40) and (3.47) as 

Qpo 



po 2 



= 2(Ve + m + l)(( + ()- 1 (4.4) 



2T° 

— = i(eve-eve + £-!)(( + c) 771 " 1 (4.5) 

The two residual Ricci and Bianchi equations now become 
££V 3 £ + UV£V 2 £ + £(2m + 1)V 2 £ - 2£V 2 £ - 2(V£) 2 - (4m + 2)V£ - 2m(m + 1) = (4.6) 
\ll 2 V ?> £ + £V 2 £((2m + 1)1+ (2m - 3)£ + UV£ + £V£) + £(V£) 2 V£ + (2m - 3)£(V£) 2 

+ (2m + 1)£V£V£+ (Am - 2)£V£- (Am 2 -Am- 3)£V£ - (Am 2 -Am- 2)^| 
- jc.c.} = 

(4.7) 

and the inequalities become 

£V£-£V£ + £-£^0 (4.8) 

£V 3 £ + W 2 £V£ + (2m + l)V 2 £ ^ (4.9) 

These two equations (4.6) and (4.7) are not independent, which can be seen as follows. 
The terms V 3 £, T> 3 £ can be eliminated from the second equation using the first giving 

\l£{plV 2 £ + (2m - l)V 2 £) + £(V£) 2 V£ + (2m - l)£(V£) 2 

+ (2m + 1)£V£V£+ (Am 2 - l)£V£ - (Am - 2)£V£ - (6m 2 - 2m - 2)£ j (4.10) 
_ j cx j = o 

25 



and it is then found that a derivative of this third equation is related to the first by 



The equations (4.6,10) and inequalities (4.8,12) for £ carry essentially the same information 
as the equations (111,112) for i given in [1], but the correspondence is not so obvious; for 
instance we note the additional explicit unknown function / given in [1] , and point out that 
in this paper there is also an additional unknown function, namely W, which is implicitly 
built into the operator V. 

There are different possible ways to exploit the redundancy noted above. From one point 
of view, it means that we have essentially only two independent real equations, e.g. (4.10) 
and the real part (4.6); from another point of view we can think of the complex equation 
(4.6) as the main equation with (4.10) as a supplementary equation, which is essentially a 
special first integral of (a part of) (4.6). We could, at this stage, simply write V = W-^ 
with v as coordinate, and attempt to decouple and simplify the two equations, in coordinate 
form, making use of the redundancy and coordinate freedom; we shall show in a subsequent 
paper how this can be done, continuing on from the equations (111,112) for £ given in [1]. 
But in this paper we shall prefer to exploit the redundancy and simplify the equations 
(4.6,10) further, yet remain within the operator formalism; we shall then show that the 
final coordinate choice and decoupling follows in a very natural mannner. 
We shall consider the pair of equations (4.6,10) as our 'master equations', and will simplify 
them in different ways in the remainder of this paper, as well as pointing out how other 
approaches from these equations will enable us to retrieve the results of earlier workers. The 
equations and inequalities now have no terms which are functions of the three coordinates 
R, C, C, being functions only of the coordinate candidate u. This of course means that when 
we introduce explicitly our fourth coordinate we will be dealing with a pair of ordinary 
differential equations in one real variable. 

4.2. The master equations rewritten as a complex chain of first order equations. 

We will now rewrite (4.6) and (4.10) in a more concise and manageable form. Defining 



V [(4.10)] = (2m - 1 + V£) [(4.6)] - (2m - 1 + Vi) [(4.6)] 



(4.11) 



In the same way as above, inequality (4.9) can be simplified to 



£V 2 £ + (VI) 2 + (2m + 1)V£ + m(m + 1)^0 



(4.12) 



A =£{V£ + 2m - 1) 
2A 

A =VX + — + (m - l)(m - 2) 



(4.13) 
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we find that the differential equation (4.6) can be written in the simple form 



VA = 2A/£ (4.14) 

while the differential equation (4.10) yields simply 

AA - (to - l)(m - 2)(2£ + A) = AA - (to - l)(m - 2)(2£ + A) (4.15) 

(We can divide by £ since a non-zero £ is guaranteed by the inequalities.) 

The inequality (4.8) remains in the same form, while the inequality (4.9) becomes very 

simply 

VA^O (4.16) 
So an equivalent presentation of these equations is, 
VA =2 A/ 1 

V£ =1 - 2m + (B + 2(m - 1)(to - 2)/)/£(A - (to - l)(m - 2)) (4.17) 
VB =AA - (to - 1)(to - 2)(A + A - to(to + 1)) 

where £(= AA - (to - 1)(to - 2) (A + 2£)) is a real function. 

(The denominator in the second equation is non-zero, since A cannot be constant, because 
of (4.16); nor can £ be zero, as we have noted above.) 

At first sight this set of first order (operator) differential equations — under the obvious 
substitutions of £ from the first into the second, and B from the resulting equation into 
the third — appears to reduce to a real differential operator equation of third order for 
the complex function A; however, a complication is that since B on the right hand side of 
the second equation is real, this equation also carries the implicit information 

£A(V£ - 1 + 2to) - (to - 1)(to - 2) (V£+l + 2m) 

- - - (4-18) 

= £A(V£ - 1 + 2to) - (to - 1)(to - 2){V£ + 1 + 2m) 

which complicates the deceptively simple structure in (4.17). Of course we also still need 
to introduce the fourth coordinate explicitly through the operator T>, which will introduce 
another function; ultimately we will need to decouple the real and imaginary parts of A. 
We note that we have exploited the redundancy (4.11) to, in effect, replace the 7 real 
equations (4.13,14,15) in the 6 real unknowns £, A, A by the 5 real equations (4.17) in the 
5 real unknowns £, £?, A; in particular, it is emphasised that the introduction of the real 
function B means that the set of five equations (4.17) has no redundancy. 
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4.3. Introducing different coordinates to achieve decoupling. 

The above set of equations reveal a structure built around the complex function A which 
we know depends on only one variable. In Section 3 we chose the coordinate candidate 
u as coordinate, and we could of course do the same for the alternative set of equations 
(4.17). In that coordinate the operator V = and substitution into (4.17) will lead to 
two real differential equations — of order three and two respectively — for the complex 
function A; although these are of a simpler structure than their counterparts obtained in 
Section 3, the difficulty of obtaining decoupling — after the separation of A into real and 
imaginary parts — still remains. 

However, we shall illustrate two approaches respectively in the next two sections, whereby 
we can obtain decoupling in a natural way. 

With our additional freedom of choice for the fourth coordinate, there is the possibility - 
suggested both by the structure of the equations (4.6,10) and the absence of u explicitly 
in these equations — to choose the fourth coordinate in terms of A. It would of course be 
very attractive if we could choose A itself as the coordinate; unfortunately A is a complex 
function of the real variable u and so A itself cannot play the role of the fourth coordinate; 
however, some real combination formed from A (e.g. K(A), arg(A), ...) would be a possible 
choice for the fourth coordinate. With such a choice the decoupling problem for complex 
A disappears, since one part of A is now the coordinate, and the other part is the real 
dependent unknown; there is also a reduction in order. (It might appear that we would 
also be left with a supplementary first order equation to solve for u as a function of A, at 
the end of the analysis. In fact, as noted above, u actually does not occur explicitly in any 
of the spin coefficients or Killing vector components and so it would not even be necessary 
to calculate u explicitly.) 

On the other hand, if we could rearrange our equations (4.17) into a closed chain of real 
equations in real functions then we could choose as our fourth coordinate the real function 
at the top of the chain. 

5. A class of real third order differential equations for all values of m. 

Let us label our fourth coordinate, v, so that we can write 

V = w4- (5-1) 

dv 
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where W{= f^) is a function of v. The second of the three equations in (4.17), with the 
first substituted, becomes 



4A 2 (A-(m - l)(m - 2)) (A + (In W), V A) 

+ A 2 (VbA + 2(m - l)(m - 2)(2to + 3)A - 2(1 + 2m)A 2 ) = 

while the third can be written as 



(5.2) 



(WB) tV -WBQnW), v =AA-(to-1)(to-2)(A + A-to(to+1)) (5.3) 

where ' denotes differentiation with respect to v. 
From (5.2) we obtain 



(lnW), v = A (3 + A 3 (3 

WB = -4(A - (to - 1)(to - 2))A 2 A/3 - 4 (A - (to - 1)(to - 2))A 2 A/3 

where 

= (2A 3 A-(m - l)(m - 2) (2A 2 A - (2m + 3)AA 2 ) - (1 + 2m)A 2 A 2 ) 

1 2AA((A - (m - l)(m - 2))A 2 A 2 - (A - (m - l)(m - 2))A 2 A 2 ) 



(5.4) 



(5.5) 



(We note that we can assume the denominator of (3 is non-zero; otherwise we obtain a flat 
space solution — with inequality (4.16) violated — or the Hauser solution [19,20,21].) 
Substituting the two equations (5.4) into (5.3) gives a real third order equation for the 
complex function A of the real coordinate v, 

(A-(m - l)(m - 2)) (4A 2 A0 - A 3 p 2 - A 3 f3j}) + 2(2m + 7)AA 2 ^) 

+ (A - (to - 1) (m - 2)) (^4A 2 A(/3 - A 3 /3 2 - A 3 (3(3) + 2(2m + 7) AA 2 (3^) ( 5 - 6 ) 
= -AA + (to - 1)(to - 2) (A + A - m(m + 1)) 

But we have still the freedom to choose v; and providing we choose v as a real function 
of A, and make an appropriate choice for the dependent variable — essentially choose it 
as a second, independent, real function of t — then we can obtain the resulting equation 
which is of third order. 
As an example, we put 

A = v + iX (5.7) 
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and write out the equation for the special case m = 1. When we substitute (5.7) into (5.5), 
for this case, we obtain 

0=(v + iX) 2 (2i{v + iX)X - 3(1 + iX) 2 ] 

/2(1 + X 2 ) [{v - iXf{l + iX) 2 - (v + iXf(l - iX) 2 } 
where ' denotes differentiation with respect to v. 

When this is in turn substituted into (5.6), which for m = 1 simplifies to 

(v+iX) 2 (l + iX) (2(v + iX) (0 - (1 + iXfp 2 - (1 - iXffifi) + 9(1 + iX)0} 

+ (v- iX) 2 (l - iX) (2(v - iX) (0 - (1 - iXf/3 2 - (1 + iXf 00) + 9(1 - iX)0} 
= -(v 2 +X 2 )/2 

(5.9) 

Clearly there is no decoupling problem, and we have a real third order equation for the 
real function X of the real coordinate v. Picking out the third order expression, we can 
write it explicitly as 

(AA - AA) = -2iX (5.10) 

We emphasise that once we have solved this single equation (5.9) — subject of course to 
the inequalities (4.8,16) — then the problem is essentially solved; once X is obtained we 
can write down A from (5.7), obtain (lnVF) iW from (5.4), and W by integration; hence we 
can obtain I from the first of (4.17) by substituting for W and A. 

So, we can then write down our table for the fourth coordinate v, 

=0 

~ / (ill ~ I — 

l>v=-i>u = w\P"\(C + cr^ ({U2) 

~ fill ~ — 

dv =—du = wp°e(( + c)" 1 

du 

which gives us, in the coordinate system R, (, £, v, the GHP operators when acting on 
zero-weighted quantities as, 

p = p(#) d R + p(C) d c + p(C) d~ c + $(v) d v = 2\P°\- 1 d R 

i>' ' ={p° • + p ,0 )\PT 1 dR + w\p \(c + cr 1 - m d v 
+ p o d c + wp°£(( + ()- 1 d v 

d' =(^(P°E^ + WP°l(C + C)" 1 ^) - RP°(Wl+m + 1)(C + C)" 1 )^ 

+ P°d- C + wP°l(c + c)- 1 d v 
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(5.13) 



where E° is given by (4.5) and £, W are obtained as described above. 
The Killing operators remain as in (3.50,51). 



Finally, we emphasise that the reduction to a real third order equation of a real function is 
not only for the particular coordinate choice v = §ft(A); we could have chosen as coordinate 
any real function of A (e.g. , 3(A), arg(A), |A|, ...) and with an appropriate choice of depen- 
dent variable, achieved an equation of the third order. So there is a whole class of real third 
order equations of a real unknown which can easily be found from (5.6). As noted above, 
the leading term in (5.6) has the form (A A — AA), and this raises the question whether 
one of these alternative coordinate choices would enable third order terms to cancel, and 
yield an equation of second order. For the special coordinate choices just mentioned it is 
easy to see that such reduction does not occur, and indeed we shall show in a subsequent 
paper that there is no possible way to choose our coordinate as a real function of i in order 
to achieve reduction of order. Of course this does not mean that there may not be some 
other choices of coordinate for which reduction of order can occur. 

6. A concise third order equation for the case m = 1. 

In this subsection we reduce the system of complex equations (4.17) to a real system, 
which form a closed chain of real first order (operator) differential equations. So as not to 
obscure the technique by details we consider the special case m = 1. Beginning with the 
one real equation from (4.17), 



VB = C 



(6.1) 



where 



B = AA 



C = AA 



(6.2) 



we use the other equations in (4.16) to obtain the two real equations 



VC = 2CE 



(6.3) 



and 



B B 1 1 



(6.4) 



VE = 



AP A£ 3 P i 2 



where 




(6.5) 
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with B, C, E ^ 0. 

We could obtain the next equation in the chain by operating with V on the right hand 
side of (6.4); but the calculations are shorter if we introduce the complex function F by 

F = £ 2 A (6.6) 

and note that not only can we write out the right hand side of (6.4) in terms of real £?, C, E 
and complex F, but we also have, from (4.17), the very simple and useful result 

VF = 2B (6.7) 

Since B is real we can put 

A + ik = F = £ 2 A (6.8) 

where A is a real function, but k is a constant. So we have now obtained the following 
closed chain of real first order equations of real functions: 
DA = 2B 

VB = C 

VC = 2CE (6.9) 



VE = -^-^ [AE ± Je 2 -aJ^-^) +E 2 -2j ° 



A 2 + k 2 \ V V A 2 + k 2 J V A 2 + k 2 

Although we now have four real differential equations compared to the (equivalent of) 
five real differential equations in (4.17) we have not lost any information; the missing 
differential equation is simply Vk = 0; and by using only real functions all the information 
is explicit within the chain. Substituting the first three equations of (6.9) into the fourth 
one will clearly give a fourth order operator equation; however, an appropriate choice of 
coordinate will reduce the corresponding coordinate equation to third order. Let us choose 
A as our fourth coordinate so that from the first equation the operator V is given by 

V = 2B^- (6.10) 
dA 

and substituting from the first three equations of (6.9) in the fourth equation gives a real 
third order differential equation for B, a real function of A, which can be presented in the 
comparatively concise form, 

2SSS + S 2 S - 3SS 2 

(6.11) 



s 
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S \ s 



S 2 S 2 -4SS 2 \ — — -r -25 



where ' denotes differentiation with respect to A and 



S = B' 



(6.12) 



The constant parameter k can be absorbed by the following relabelling 

A =A/k 
S=S/k 

giving 



(6.13) 



2SSS +S 2 S -3SS 2 



s 



(A 2 + 



— (AS§+^ S 2 S 2 -4SS 2 < 



s 



A 2 + l 



-2S 



S 



A 2 + l 
(6.14) 

where ' now denotes differentiation with respect to the new coordinate A. 

It is clear that once S is obtained the problem is essentially solved: C, E can be obtained 

from (6.9) by differentiation, and can be combined algebraically to give £, 



2Vk 



ss 



+ 



i 



ss 



s 



A 2 + l 



(6.15) 



The table of GHP operators in this coordinate system, R, £, £, A is obtained from (5.10) 
where v is replaced by A and then by A, and W is replaced by B, and then by S. 
We emphasise that the inequalities (4.8,16) must also be satisfied; this constraint prevents 
us from making the simple choice k = 0, because it violates the first of these inequalities. 

We could have chosen other coordinates e.g. B or C or combinations, and in the same 
manner obtained alternative third order equations. 

Finally, we note the relative simplicity of the version of the final equation obtained here 
as compared to the version, for m = 1, in the last section. 



7. Summary. 

The original purpose of this paper was not primarily to try and make significant new 
progress in the NT problem; it was rather to set out in detail new insights into the GHP 
operator-integration approach — whose principles are applicable in very general contexts 
- and then to illustrate these insights by applying them to the NT problem. But, as well as 
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the NT problem providing an ideal laboratory for the demonstration of our GHP operator- 
integration approach, we have been able to obtain new insights into the NT problem itself 
— extending existing results and suggesting new approaches. 

The major new insight regarding the GHP operator-integration approach is the important 
fact that, although spacetimes admitting Killing vectors do not immediately supply — from 
the spin coefficients and Riemann tensor components — the four functionally independent 
scalars which integration within the GHP formalism demands, the missing scalars may 
be supplied in a very simple and natural way — by also using the (tetrad components 
of the) Killing vectors. The relationship in the GHP formalism between Killing vectors 
and functionally independent scalars, as well as that between homothetic Killing vectors 
and separability, are very important topics which we have only touched on here; deeper 
implications for the GHP formalism will be developed in a separate paper. 
An additional insight regarding the GHP operator-integration approach is that, even after 
we have reduced the problem to the residual ordinary differential operator equations, there 
are advantages to remaining within an operator formalism rather than immediately using 
coordinates explicitly. We have shown that, by constructing a closed chain of real first 
order ordinary differential operator equations, we can avoid the difficulties associated with 
decoupling, and can make the final coordinate choice in an efficient manner. 

We chose the G2 case of the NT problem as our illustration because of its non-trivial nature 
and because of the wide range of mathematical procedures necessary for its simplification. 
The intention was that the advantages and strengths of our method be thoroughly illus- 
trated and tested. Much work has already gone into this particular problem, and many 
insights have already been gained — although often in a rather narrow manner specific 
to a particular formalism and coordinate system. In this paper, we have been able to 
understand the procedures of separability, redundancy, decoupling and reduction of order 
in a very general manner, in the context of the GHP formalism; we have also seen how 
it is possible to retain our coordinate choice to the very last step when it can be used to 
simplify the final equations in the most advantageous way. This we believe is the strength 
of this GHP operator approach; separability, redundancy even decoupling can be exploited, 
still retaining some coordinate freedom, and then at the very last stage when everything 
has been reduced to the decoupled residual ordinary differential operator equations this 
freedom can be exploited in such a manner as to present the final equations in the most 
reduced and/or manageable form. 
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Other workers, [15,16,17,18] who have reduced this problem to third order real differential 
equations have used very different formalisms, and special techniques, and there is no 
simple relationship between them. However, it is possible to obtain these very different 
equations from our master equations: 

(i) The form of the residual real third order ordinary differential equation for h(ip) obtained 
by Ludwig and Yu [17] follows in a very concise manner from our master equations (4.6,10), 
by making the special choice of coordinate cp = arg(£) and of independent variable h = 

We note, that unlike in the approach in Section 5, the particular coordinate choice 
is crucial; most other coordinate choices built on i give fourth order equations. 

(ii) For arbitrary values of the parameter m, Herlt [16] has first of all reduced the problem 
essentially to a real third order ordinary differential equation for a complex function g, 
which he then transforms to a first order integro-differential equation for the real function 
<p(u) by defining the coordinate u by 



Finally, he is able to obtain a real third order ordinary differential equation for a real 
function x(f) defined by 



It is easy to see that Herlt's function g corresponds essentially to our function A, and 
so we can retrieve his result by first rearranging our master equations (4.6,10), and then 
following his procedure. However, it is also possible to obtain other third order equations 
by constructing different coordinate choices around A (equivalently Herlt's g) than the 
choice (7.1). 

For the special case m = 1 (corresponding to N=2, in the notation of [16]) Herlt has 
obtained a particularly concise form for the residual real third order differential equation, 
just as we have also, for this special case — in Section 6. Although there are obvious 
similarities in structure, some manipulation still needs to be carried out to show the direct 
equivalence of the two equations. 

(iii) Chinea [15] — using an elegant approach with matrix valued differential forms — has 
reduced the NT problem with two commuting Killing vectors (i.e. m = 0) to a single 
complex second order differential equation for a complex function. If we choose B as our 
coordinate, then after substituting the first and third into the second of equations (4.17) we 



9 = 9, u e 



.i<p(u) 



(7.1) 




(7.2) 
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obtain a complex second order differential equation for arbitrary values of the parameter 
to; for the special case to = this equation is 



(AA - 2A - 2A)(2AA - A 2 ) = 2A(2AA - 2AA + A 2 A) - B(AA -2k- 2A) 2 A 3 /2A 2 (7.3) 



where ' means differentiation with respect to B. This equation is similar in structure to, 
but is not yet as neat as that of Chinea [15]. However, for non-commuting Killing vectors 
the two special cases to = 1, 2 have the remarkably simple forms, given respectively by 



Chinea was able to transform his complex second order equation, by a number of coordinate 
changes, to a very complicated real third order ordinary differential equation for a real 
function of a real variable. We could do the same for the equations above (and also for 
the equation for arbitrary to), but the resultant equations seem more complicated than 
the versions obtained in the previous sections. 

(iv) Finley et al. [18], beginning from the point of view of groups of point transforma- 
tions, also reduce the problem to one real third order ordinary differential equation for a 
real function of a real variable q(y). They also define a function which is second 

order in q and point out that the complicated third order equation can be rewritten as a 
comparatively simple, but interesting, first order expression for the real function *&(y) - 
spoiled only by one additional term explicit in q(y). Mcintosh [23] has also found a similar 
interesting first order expression, again, unfortunately with one additional complicated 
term. As Finley et al. point out, these first order expressions have suggestive symmetries, 
which may lead to further significant simplifications. Essentially the idea being followed 
is to write the very complicated third order equation as a coupled pair of equations; in 
particular one of these equations should be simple enough to enable solutions for it to be 
found, yet complicated enough that the solutions are not trivial, and also that the second 
equation in the coupled pair should be significantly simplified. 

In fact, our versions of the master equations in both (4.17) and (6.9) are in a form which en- 
ables us to experiment in precisely this manner. For instance, let us choose A as coordinate 
and rewrite the last equation in the set (6.9) as 



AA + AA = (2A - 5AA)AA /4A 2 



(7.4) 



AA + AA = (6A - Bkk)kk /4A 2 



(7.5) 



2X, A -X 2 + 



AX 



= Y 



(7.6) 



A 2 + k 2 
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where we have relabelled X = E/B, and where Y is given from (6.9) by 



Y = 



XC 



± 



k 



X 2 



4 / C 



B 2 y a 2 + k 2 



2 / C 



(7.7) 



2B 2 



A 2 + k 2 



B 2 \ A 2 + k 2 



From this equation we can easily find an expression for X in terms of B, C, Y and the coor- 
dinate A. If we differentiate this expression twice we get two new expressions respectively 
in Y^a, Y, £?, C, and Y^aa, Y,a, Y, B, C; by eliminating £?, C between these three expressions 
we are left with an expression for X in terms of Y, Y^, Y^a and the coordinate A i.e. an 
expression for X which is second order in Y . 

The structure of the left hand side of (7.6) is very similar to the interesting first order 
equations obtained respectively by Finley et al. and Mcintosh. 

(v) Finally we point out that the general third order equation which we obtained in Section 
5 can also be obtained using the coordinate method of our earlier paper, [1] . By a suitable 
coordinate transformation, and appropriate combination of equations (111) in [1], we can 
obtain [5.6], after a lengthy calculation. 

We shall show in a further paper the full details of how all of the existing equations, 
together with some simpler new ones, can be obtained from our master equations. 
It is clear that there are many versions of a final third order equation, and so it is certainly 
possible that more manageable ones exist; there is even the possibility that a second order 
equation may still be found. We believe that the approach in this paper, which has the 
power to give us a unifying picture of such different approaches, supplies us with very 
powerful equipment for this ongoing task. 

APPENDIX: SIMPLIFICATION OF GHP EQUATIONS IN INTRINSIC 
TETRADS. 

A.l. Spacetimes admitting a Killing Vector. 

We emphasise that in [12] the system of equations given was for the presence of a homo- 
thetic Killing vector, while in [8,9] only the presence of a Killing vector is considered. 
In this subsection we are only considering Killing vectors, so when we refer to the equa- 
tions in [12] we understand that we have made the substitution = which reduces the 
conformal Killing equations to the Killing equations. 

In [8,9] the Killing equations consist of the seven complex and the two real equations 
(A.8a,b,c,9a,b,10a,b,13,14) while in [12] the Killing equations consist of the eight complex 
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and four real equations ((21), (21)', (21)*, (21)'*, (22), (22)', (22)*, (22)'*, (23), (23)', 
(23)*, (23)'*). 

The discrepency in the number of equations is easily explained when we realise that (A13) 
in [8,9] has simply been split in [12] into the two equations (22)' and (22) by the introduction 
of an arbitrary real quantity V\ similarly (A14) in [8,9] has simply been split in [12] into 
the two equations (22)* and (22)'* by the introduction of an arbitrary real quantity V*. 
Finally in [12], (21)'* is simply the complex conjugate of (21)*; in [8,9] only (A. 8b), the 
counterpart to (21)* in [12], is given explicitly. 

But there is a second apparent discrepency. In [12] there occurs in (23) and (23)'*((23) / 
and (23)*) an arbitrary complex quantity Q (Q f ); no such arbitrary quantities occur in 
the counterpart equations (A9a,b) and (AlOa,b) in [8,9]. However, it is easy to see that 
when we apply the vacuum Type N restrictions with the particular tetrad gauge used in 
[12] to the complete system of equations in [12], that it follows that 

Q = 0=Q' (All) 

(In fact, in [12] this calculation was carried out for the more general system — admitting 
a homothetic Killing vector — resulting in (AI.l). 

Therefore, the apparent discrepency between the two systems of equations in [12] and [8,9] 
respectively, in the presence of a Killing vector is resolved — for the particular class of 
spacetimes under consideration in this paper. 

However, we know from Held's argument that such simplifications can be made in much 
more general circumstances. This is a very fundamental property for spacetimes admitting 
homothetic Killing vectors, and has very important consequences for the GHP formalism. 
We will deal with this topic elsewhere; here we will just illustrate other more general cases 
where this simplification can be made. 

(i) For instance, suppose we choose l % along the principal null direction of the Weyl tensor, 
so that 

^o = 

Then we can deduce that 

Q = 

provided that at least one of the ^l/n's (n = 0, ...,4) is nonzero, by an iterative argument 
using in order (43), (44), (45), (44)', from [12]. (See [22] for a related type of discussion.). 



38 



Similarly we can choose n l along the principal null direction of the Weyl tensor, so that 



#4 = 



and this is compatible with 



Q' = 



(We note that the choices #3 = = #1 also lead to the same result.) 
(ii) For algebraically special spaces if we choose l l such that 



W = = #1 



then, as in (i), 



Q = 



We can use the null rotation about l l to put 



Q' = 



In an analogous manner we could fix l l with respect to the Ricci tensor etc., [22]. 
A. 2. Spacetimes admitting a Homothetic Killing Vector. 

We now consider the set of equations in [12] in the presence of a homothetic Killing vector, 
and as noted above, it has been shown in [12] that (AI.l) follows when we apply the 
vacuum Type N restrictions with the particular tetrad gauge used in [12]. Therefore, for 
this special case, we can combine the 12 equations in [12] into the following 9 equations, 



f> a 

pda 

pdb 
P c 







— k'c — k'c 



-pc 



(p - P)c 



R'b + (p' - p')c 



(AI.2) 








-p'c 



~ 1 



pb c + pbc 



<f> - (p' + p')b - (p + p)a 
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These equations are just the homothetic Killing vector generalisations of the analogous 
Killing vector equations (3.1)-(3.9) in [8], specialised to vacuum Type N spacetimes. 
We could have derived these equations for much more general classes of spacetimes, using 
a generalisation of the geometric argument used by Held in [8]; however, this would involve 
us in deeper questions than we wish to consider here, and the above equations are sufficient 
for our purposes. 
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